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Abstract 



We examine deviations from Boltzmann-Gibbs statistics for partially equi- 
librated systems of finite size. We find that such systems are characterized by 
the Levy distribution whose non-extensivity parameter is related to the num- 
ber of internally equilibrated subsystems and to correlations among them. This 
concept is applied to Quark-Gluon Plasma formation. 

PACS numbers: 05.70.Ln, 05.40.Fb, 25.75.-q 

1 Introduction 

Non-extensive statistics introduced by Tsallis [p]] has become a topic of substantial in- 
terest and investigation in many areas of physics including condensed matter, nuclear, 
and high energy physics. The centerpiece of such considerations is the power-like Levy 
distribution 



where 1 < q < 2 and C q is the corresponding normalization constant. In the limit 
q — > 1 this reduces to the standard Boltzmann-Gibbs exponential factor, 
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In this case the parameter Ao = l//?o can be interpreted as the temperature of the 
system in statistical equilibrium. (We set the Boltzmann constant equal to unity.) On 
the other hand the meaning of Eq. ([I]) is still not perfectly clear for q ^ 1. In a recent 
article Wilk and Wlodarczyk 0] have shown that the Levy distribution naturally 
arises as the mean of the Boltzmann-Gibbs factor over a Gamma distribution of 
the temperature parameter, 1/A = 1/T, with mean 1/Ao- In the same article the 
authors have shown that the parameter q is determined by the diffusion coefficient 
and the dumping constant in the Fokker-Planck equation, i.e, the Langevin equation 
describing gaussian white noise. We wish to further elaborate on this observation by 
finding other possible origins of this Gamma distribution. 

Tsallis' statistics, initially introduced to describe fractal properties of various sys- 
tems, maintains the general structure of thermodynamics. Formal analogies between 
such statistics and the so-called g-oscillators have been established by means of q- 
calculus and use of Jakson derivatives This statistics has found applications in 
many phenomena including modifications of solar neutrino fluxes ||, fluctuations and 
correlations in high-energy nuclear collisions [|], the long-flying component of cos- 
mic rays ||, chaotic transport in laminar fluid flow 0, subrecoil laser cooling || and 
others. It has also been incorporated into fractional quantum statistics ||. 

Gross and Votyakov [[HJ have investigated the statistical properties of small sys- 
tems (systems whose size is comparable with the range of internal interactions) and 
have concluded that, in this case, the grand-canonical distribution is not appropri- 
ate and the thermodynamic limit cannot be taken. Earlier Prosper [IT] had already 
studied temperature fluctuations in finite-size heat baths and shown that the energy 
distribution of a thermometer is, in fact, binomial. It appears that finiteness of the 
system size and correlations between its parts due to long range interactions can lead 
to substantial deviations from the widely-used Boltzmann-Gibbs distribution. An 
exhaustive coverage of all work related to Tsallis' statistics can be found in Ref. |T2 |. 



In this article we examine these concepts in the case of partially equilibrated systems. 



2 Derivation of the Levy Distribution 

First we would like to see how the Levy distribution can arise for finite-size partially 
equilibrated systems and to investigate the meaning of its parameter q. We consider a 
Gibbs ensemble of finite-size system replicas each subject to the constraint of fixed to- 
tal energy, E [13] . Each system, member of the ensemble, is partially equilibrated |14| . 



It consists of finite-size subsystems each of which is in internal statistical equilibrium 
corresponding to some temperature, T. Since the size of the subsystems and of the 
whole system is finite the number of subsystems is also finite. For any member of the 
ensemble there is a number k of subsystems at temperature T (more precisely in the 
interval between T and T + dT). The ensemble mean of the number of subsystems 
at temperature T is (k). We define To, the ensemble mean of the temperature of all 
subsystems first averaged over the number of subsystems of each ensemble member. 
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This is the same as the temperature average over all the subsystems of all ensemble 
members. The more ensemble members have average temperatures close to T or 
the more subsystems over the ensemble have temperatures close to To the closer to 
equilibrium the physical system will be. We assume that the energy of the whole 
system is proportional to the mean temperature, 

E = cT 0} (3) 

where c is a constant. At the same time the ensemble mean energy of all subsystems at 
temperature T will be taken equal to the ensemble mean of the number of subsystems 
at this temperature times the energy of each subsystem, i.e., 

(E(T))=c(k)T. (4) 

Under the strong assumption that E = (E(T)) we can easily see that 

(*) = p = Ao/8, (5) 

where Ao = T and f3 = 1/T = 1/A. This condition correlates the subsystems and may 
be the result of long-range interactions within the system. The case T = T is rather 
special since, then, Eqs. (|],[5]) imply that (k) = 1 and, in this case, there can be only 
one subsystem at temperature equal to the average one. The probability distribution 
for encountering exactly k subsystems at temperature T inside one ensemble member 
can be considered to be Poissonian for those values of k reached within one ensemble 
member, 

Consequently, the probability for encountering a subsystems at temperature T within 
the ensemble is given by the Gamma distribution 

r(ar) 

with first moment a/Xo- Upon rescaling A by dividing it by a and multiplying Eq. (|7|) 
by a we obtain the modified Gamma distribution 

'-(i 

with first moment 1/Ao- This is the same function as the one obtained by Wilk and 
Wlodarczyk || starting from the Levy distribution. We may define a parameter q 
such that 
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Clearly when q — > 2 we have only one subsystem at temperature To in the ensemble 
(a = 1). The system is away from equilibrium. We note that if the first moment of 
the Levy distribution is to be finite then q < 3/2 so that complete lack of statistical 
equilibration cannot be achieved. When q — > 1, a tends to infinity; there is an infinite 
number of subsystems in the ensemble at To. In this case the system is equilibrated. 
To derive the Levy distribution we follow the inverse steps of Ref. ||. First we define 
a variable £ by means of the equation 

€ = ^- do) 

If each subsystem is attributed a Boltzmann-Gibbs factor in a variable x, exp(— x/X), 
then the expectation value of this factor over the distribution given by Eq. (§) is 



i 



/ r _1 exp 
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d£. (11) 



Even though for a particular ensemble member 1/A cannot cover the entire range 
[0, oo] this is cetrainly possible over the ensemble. We may make the substitution 
t = £ [1 + x/ (ctAo)]. This yields 

<«^>=i^( i+ ^)T^M i+ ^r « i2) 

and which, by means of Eq. @, gives the unnormalized equivalent of Eq. ([!]). In 
this derivation the finiteness of the system and subsystem sizes has been important 
in order to assume a (discrete) Poissonian distribution for k. 



3 Entropy Non-Extensivity 



We, now, consider the entropy of the system. Each subsystem being in internal 
equilibrium has a Boltzmann-Gibbs-Shannon entropy, 



S; 



Wi 



(13) 



where Wi is the number of microstates consistent with a given temperature of the 
subsystem and Pij is the probability for the jth microstate of the subsystem of tem- 
perature Tj. Since each subsystem is characterized by a g-value, q , very close to 
unity we can apply the approximation pfj -1 ~ 1 + (qo — 1) hip^. In combination with 
Eq. QTJ) this yields 
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This power-law entropy was first introduced by Tsallis |15| and upon optimization 
under the constraints of normalizability and given second moment produces the Levy 
distribution. Even though this equation was derived as an approximation its applica- 
bility is more general than that of the Boltzmann-Gibbs-Shannon entropy ||12|| . This 
form of the entropy is, generally, non-extensive as can be seen by calculating the 
combined entropy of two stochastically independent (sub) systems, 

1 



qo - 1 



1=1 k=i 



S T + S v - (q - 1) S T (15) 



where the combined system is assigned the same g-value as its components. The last 
equation can be proven by considering the product of the two individual entropies. 
Clearly as q — > 1 extensivity is recovered. If all subsystems are equally close to inter- 
nal equilibrium and are stochastically independent the entropy of the whole system 
is 

s = E^-(g -i)EE^-(g -i) 2 E E E SiSjS*-... (ie) 

i i j>i i j>i k>j>i 

A microstate of a subsystem is part of a microstate of the system. If a particular 
subsystem exists in one ensemble member but does not exist in another then the 
contribution to the probability of the subsystem microstate is zero for the latter 
member. Therefore, the above summations extend over all the subsystems that may 
appear in any ensemble member. We observe that in the case of stochastically in- 
dependent subsystems the total entropy is smaller than or equal to the sum of the 
individual entropies, the equality always holding when the subsystems are internally 
equilibrated. However, Eq. (|5|) induces correlations among subsystems which remove 
the stochastic independence. It is due to this effect that the total entropy of the 
system is non-extensive even though its parts are in (or close to) internal equilibrium. 
This conclusion extends the one reached by Landau and Lifshitz [14]] on partially 
equilibrated systems. 

In the presence of long-range correlations and if the stochastic dependence is 
relatively weak we may approximate the sum of the g-power products of the microstate 
probabilities of two (sub) systems by 



Wr Wu Wtv 

1=1 k=l p=l 



EE^M^E^ + c- (17) 



where W TV is the number of microstates of the combined system, p p the probability 
for one such microstate, C TU is a function of go that incorporates the correlations 
and g is a parameter defined below. We express the functions C TU in terms of the 
individual entropy product as 

C TV = k S T S U) (18) 
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where k = K,(qo). Using the last two equations we find that the total entropy is 

S TV = S T + S v - [g - 1 + k/ (g - 1)] S T S u = S T + S u - (g - 1) S T S v , (19) 

with g = g + n/(q — 1). When k^O the total entropy is not extensive even though 
the individual subsystems may be close to internal equilibrium. We note that k must 
be an increasing function of go — 1 to prevent the paradoxical situation of g becoming 
larger the closer g is to unity. 



4 Perspectives for Heavy Ion Physics 

We have seen that the condition of partial statistical equilibration in a finite-size 
system leads to a Levy distribution that generalizes the Boltzmann-Gibbs statistics. 
Long-range correlations among the internally equilibrated subsystems are important 
and give rise to the non-extensivity of the entropy. The parameter q that characterizes 
this situation is related to the correlations and to the number of internally equilibrated 
subsystems. It is certainly interesting that the very complicated dynamics of such 
systems could be effectively described by a single-parameter statistical distribution. 
One possible application of these ideas is in the formation of Quark-Gluon Plasma 
(QGP) in relativistic heavy ion collisions. If such an unconfined state is actually 
produced at high temperatures but relatively low densities then the fragmentation 
time of the quarks and gluons to hadrons may be comparable to the relaxation time 
of the whole system. In this case the plasma is at most only partially equilibrated. 
Thermal photon emission from it will correspond to a distribution of temperatures 
and will effectively be described by Tsallis' statistics that corresponds to the average 
temperature of the system and a q- value that exceeds 1. It is particularly interesting 
that as Walton and Rafelski have shown [16|] diffusion of the charmed quarks in 



a plasma in equilibrium can be described well only by the non-extensive Tsallis' 
statistics. On the other hand if the plasma is not formed but photons are emitted by 
a perhaps partially equilibrated hadron system the resulting enhansement of photon 
production (for q > 1) as the photon transverse momentum increases may mimic 
photon production from a QGP. 
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